Proof of the thermodynamical stability of the E' center in SiC>2 
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The E' center is a paradigmatic radiation-induced defect in SiC>2 whose peculiar EPR and hy- 
perfine activity has been known since over 40 years. This center has been traditionally identified 
with a distorted, positively-charged oxygen vacancy Vq. However, no direct proof of the stability 
of this defect has ever been provided, so that its identification is still strongly incomplete. Here 
we prove directly that distorted Vq is metastable and that it satisfies the key requirements for its 
identification as E', such as thermal and optical response, and activation-deactivation mechanisms. 

PACS: 61.72.Bb, 71.55.Ht, 61.72.Ji 



Understanding defects in solids is a key factor in im- 
proving device performance and materials quality. Defect 
identification combines experimental observation and 
theoretical predictions, and a major ingredient in this 
process is the theoretical justification of the thermody- 
namical stability of the defect. If this information is miss- 
ing, the identification is arguably incomplete or uncer- 
tain. Surprisingly, this is the case for the E' center [jj],D , 
a radiation-induced defect first observed experimentally 
as early as 45 years ago 0] in Si02 , a material of prime 
current importance in microelectronics and fiber optics 
. E' is traditionally identified with a positively-charged 
distorted oxygen vacancy Vq H , with important support 
from calculations of hyperfine couplings [pi and optical 
activity 0; its thermodynamical stability, however, was 
never theoretically proven, and the mechanisms involved 
in its activation and deactivation are still debated. Here, 
using the ab initio theory of defect formation in solids J7| , 
we show that the conditions for E' stability are naturally 
realized in stoichiometric or neutron- irradiated Si02 , and 
conclusively put on firm ground the identification of the 
E' center. 

Our model of the stability of E' is based on two native 
defects: the oxygen vacancy Vo and the oxygen inter- 
stitial Ov The motivation is that E' is observed chiefly 
(though not only) in neutron-irradiated material jj]-[|, 
where Vo and O; are both abundant. Indeed, we find 
that it is their concurrent presence that produces the 
conditions for the existence of E', in neutron- irradiated 
as well as non- neutron-irradiated material. In the former, 
vacancies Vo and interstitials Oi are essentially produced 
in pairs by knock-on, kick-out events; in the latter, they 
form in thermal equilibrium and, as it turns out, in sim- 
ilar concentrations. As will become apparent below, our 
argument on E' applies to both cases without modifica- 
tions, except for the fact that neutron-produced defects 
appear in concentrations determined by the irradiation 
dose, whereas the concentration of thermally formed de- 
fects depends on formation energies, which can be di- 
rectly predicted. 

At a growth temperature T g and with N s available 



sites, the equilibrium concentration of a defect is D= 
N s exp(— Ff orm /kBTg). The formation free energy Ff orm 
= Ef orm - T Sf orm depends on the chemical poten- 
tials of atoms added or removed, on the defect charge 
state, i.e. the charge released to or captured from the 
thermodynamic reservor constituted by the surrounding 
crystal, and on the electronic chemical potential fi e of the 
latter. Once the formation energies of the all relevant de- 
fects (vacancy and interstitial in our case) are known, the 
defect concentrations and the chemical potential fj, e are 
determined self-consistently, subject to charge neutral- 
ity, as detailed in H. A specific defect configuration or 
charge state is predicted to exist if its formation energy is 
lower than that of all other defect states for some value of 
/i e . Also, the defect is metastable if a non-zero energetic 
barrier prevents its deactivation or disappearance into 
other lower-energy configurations of the same defect, or 
recombination with other defects. The formation energy 
for our defects in charge state Q reads 



Eform(Q) = E t ° t (Q) 






add' 



Q[i e + M{Q) + P (1) 



where EJ^ t and E™ t are the total energies of the de- 
fected and undefected system, respectively, fi e is the 
electron chemical potential (equaling the Fermi level Ep 
in our T=0 calculations), M(Q) is the defect-dependent 
multipole correction for charge state Q of Rcf. ||, P = 
Ho for Vo and P = /io for Oi, and no is the oxygen 
chemical potential. The latter is fixed to stoichiomet- 
ric conditions, i.e. at the center of its variation range 
Mmoi/2 + AH/2 < /io < Mmoi/2 determined by the to- 
tal energy /i mo i of the O2 molecule, and the calculated 
formation enthalpy AH of Si02- Ionization energies, i.e. 
the energy needed to promote (say) an electron from the 
valence band to an empty acceptor level, are defined via 
total energy differences of different charge states. Forma- 
tion entropies are beyond the scope of the method used 
here; plausible estimates are used when needed. 

Energies and forces are accurately calculated from 
first-principles within density-functional theory in the lo- 
cal approximation, using the ultrasoft pseudopotential 
plane- wave method as implemented in VASP H . An iso- 



lated defect is simulated in periodic boundary conditions 
via the repeated supercell approach. We use crystalline 
a-quartz SiC>2 supercells of tetragonal symmetry, com- 
prising 71 to 73 atoms and of linear dimensions 18.49, 
16.02, and 20.44 atomic units (theoretical lattice param- 
eters Q, matching experiment to about 0.5%). Atomic 
geometries of the defects are optimized for each Q (ob- 
tained by removing or adding electrons as appropriate) 
until all residual force components in the system are be- 
low 0.02 eV/A. No symmetry restriction is imposed on 
geometry optimization. Improving slightly on the setting 
of Ref. ||], a (222) mesh is used Q for k-space summation 
(4 special points in the supercell Brillouin zone) . Our use 
of a crystalline SiC>2-based model of the E' defect, which 
is observed both in amorphous and crystalline phases, 
is justified by the closely similar behavior of several E' 
variants in crystalline and amorphous Si02 in experiment 
|p|,m as weu as theory ||||. In addition, the simulated 
structure of amorphous silica |l2[ deviates moderately 
from that of crystalline a-quartz Si02- 
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FIG. 1. Minimum formation energies (eV) of the oxy- 
gen vacancy (dash-dotted line) and split-interstitial oxygen 
(dashed) as a function of the Fermi level, and calculated Fermi 
levels position (vertical solid). Fermi level zero is the valence 
band top. 

Figjj] shows the ground state formation energies of the 
interstitial O^ and the vacancy Vo- The vacancy acts 
as a double acceptor (double donor) in extreme n (p) 
doping conditions, but it is in fact neutral for most of 
the Fermi level range. As to the interstitial, we find 
that an oxygen atom initially placed near the center of 
the hexagonal channels of quartz, relaxes sideways to- 
wards the helical chains and, after overcoming a small 
(~ 0.2 eV) barrier, it stabilizes into a split-interstitial 
(s-Oj) configuration with the nearest bridging oxygen in 



the helical tetrahedra chain, with ~ 1 eV of energy gain 
with respect to the starting site. Details will be dis- 
cussed elsewhere: here we note the main consequences of 
this result: a) the stabilization in the s-Oi configuration 
(which is found to be quite close to that suggested in Ref. 
]13[) prevents Oj from easily recombining with vacancies, 
because the detachment from the split-interstitial config- 
uration costs about 1.2 eV; b) as seen in Fig.0, s-O,; is 
a (negative-I/) deep double acceptor with second ioniza- 
tion energy at 3.3 eV above the valence band top Ev. 
The Fermi level, calculated as in M, is pinned at Ep = 
E.„ + 3.3 eV (vertical solid line in Figjl]). Thus a con- 
sequence of vacancy-interstitial pair formation, is that 
moderately p-type conditions are achieved; this is indeed 
often observed in irradiated samples p|,|ll|. In the ab- 
sence of s-Oi, the Fermi level would be at midgap, E„ + 
4.4 eV. From the formation energies one can estimate the 
chemical concentrations of s-0 2 ; and Vo: for a typical T g 
of 1500 K, and assuming a reasonable formation entropy 
Sform= 5 ke, the concentration of both defects is ^10 14 
cm" 3 . This figure matches well E' concentrations mea- 
sured JO] in non-neutron irradiated samples after UV, 7, 
or X illumination fall, and therefore supports the hypoth- 
esis that the vacancy is the parent defect of E' . Neutron 
irradiation of course produces dose-dependent pi, typi- 
cally much higher concentrations (> 10 19 cm" 3 ). The 
key point is that, since the concentrations of vacancies 
and interstitial are essentially the same in both cases, 
and because only the interstitial has electrical activity 
(via its double acceptor level), the Fermi level is pinned 
at the same value in both cases, so that the theory of 
E' stability and activation discussed below applies iden- 
tically to both cases. 

It appears from FigJl] that the +1 charge state of the 
vacancy, Vq, is not among the thermodynamically sta- 
ble ground states of the defect. Therefore, if this state 
of Vo is to be identified with E', it must at least be 
proven metastable; if it is metastable, a mechanism for 
its creation starting from the ground state (the neutral 
vacancy) must be identified. As to the first point, since 
the +1 vacancy has a formation energy that increases 
linearly with the Fermi level while the neutral one re- 
mains constant, Vq or a distorted variant thereof, may 
only be metastable in a limited range of that variable: we 
show below that the Fermi level pinning due to split-O; 
produces naturally the conditions for the metastable exis- 
tence of E'. As to the second point, experiments indicate 
that E' is activated by ionizing radiation pp,pd| such as 
7, X, or UV photons shone onto vacancy-containing sam- 
ples, or concurrently with neutron irradiation (causing 
knock-on vacancy creation). Indeed, since the vacancy 
ground state is the neutral undistorted configuration, the 
distorted +1 state (alias E') can only be accessed by ion- 
ization of V Q : our picture provides consistently such ac- 
tivation mechanism. 

We proceed to study the behavior of the vacancy when 



subjected to the undistorted-to-puckered transition as 
proposed in earlier studies |Q pj. In accord with the re- 
sults discussed above, we fix the Fermi level at E„ + 
3.3 eV. The creation of a vacancy starting from the per- 
fect lattice results in moderate local distortions in both 
the neutral and +1 charge states. The puckered config- 
uration is obtained by moving one of the two vacancy- 
adjacent Sii and Si2 atoms (specifically the "long-bond 
HD ^2) away from the vacant site, and pushing it across 
the basal plane of the incomplete tetrahedron centered 
on Si2 itself. When Si2 pokes through this triangular 
constriction, it gets strongly and suddenly bound to a 
backbonding oxygen, 0&. Upon completion of the distor- 
tion, Si2 regains 4-fold coordination, and backbonding Of, 
becomes 3-fold coordinated (see also |§-|§]), while it was 
originally 2-fold coordinated as all tetrahedron-bridging 
oxygens in Si02- Sii remains instead 3-fold coordinated: 
in the +1 charge state, its dangling bond is half- filled, 
and causes the observed jy.2] and predicted || EPR sig- "~" 
nature which identifies E' [ 2 1 . 

It may appear at first sight that the puckering distor- 
tion should be symmetric in Sii and Si2 ■ This is not the 
case, however, because of the intrinsic asymmetry of the 
quartz structure. As already noted earlier on H], only 
Si2 finds the backbonding oxygen 0& in the correct posi- 
tion. This applies largely also to amorphous silica, whose 
structure is moderately different at the local level from 
that of quartz |l2[ . The backbonding oxygen is therefore 
the main stabilizing agent of the E' defect. 

The total energy of the system in charge state Q is cal- 
culated as a function of the separation between Sii and 
Si2- Only the modulus dsi-Si of the Sii-Si2 connecting 
vector is constrained, and all other degrees of freedom are 
fully relaxed: the minimum energy path is thus mapped 
out for the undistorted-puckered transformation in the 
constrained-dsi-si configurational subspace. In Fig. we 
display the full level diagram for the neutral, +1, and +2 
charge states of the vacancy as a function of the puck- 
ering distortion, quantified by dsi_Si- All energy curves 
depend on the Fermi level through Eq.l; they can be 
directly compared on the same energy scale because Ep 
is fixed at the natural value of E„ + 3.3 eV determined 
above. 

The outstanding feature of Fig. is that for the natural 
Fermi level of stoichiometric or neutron-irradiated silica, 
the candidate E', i.e. distorted Vq, is indeed the stable 
defect state for the distorted geometry. We stress that 
the Fermi level position is essential here: if Ep were at 
midgap, the +1 curve would be 1.1 eV above its position 
in Fig. . Then E' would be unstable towards magnetically- 



inactive Vq- Globally, E' is metastable with a confining 
barrier of 0.8 eV. The barrier to enter the metastable 
state is 1.1 eV, and the undistorted state is lower than 
the distorted by 0.3 eV [this difference to Ref. H, where 
the distorted state was lower by the same amount, is 
possibly due to our improved k-sampling]. Clearly, in 



the absence of excitation, Vq will remain trapped in the 
metastable E' state and will show EPR activity. When 
thermally activated to overcome the barrier, the puckered 
center will transform into undistorted, and by electron 
capture it will become neutral. Therewith, E' disappears 
permanently, and so does its magnetic activity, because 
Sii and Si2 combine their dangling bonds to bind into 
a dimer P~p|,^0[. (The same deactivation route is not 
readily available for the level ordering of Ref. M , which 
implies that i) E' remains activated at equilibrium, and 
that ii) concurrent barrier jump and electron capture are 
required to quench it.) 
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FIG. 2. Total energy of neutral, positive (E'), and dou- 
bly positive oxygen vacancies as a function of the puckering 
distortion, at the calculated Fermi level value. 

The calculated puckered-to-undistorted barrier is now 
compared with an estimate extracted from the measured 
relative drop in E' population upon isocronal thermal an- 
nealing pd[ in irradiated silica. In the simplest model, 
the distorted-state population No diminishes by a factor 
of N/No = exp (-R(T) r) upon annealing over a time r at 
temperature T, with escape rate R = v exp (-AF/keT), 
with AF the free energy barrier for escape from the puck- 
ered state, and v an average vibrational frequency in 
that state. Using the data |1J] for E' and assuming 
v ~ 50 THz, we obtain an experimental AF^l.l eV, in 
reasonable agreement with our calculated AF^ AE=0.8 
eV; account for the transition entropy, which we neglect, 
should further improve agreement since the transition 
state (through the tetrahedron basal plane) is severely 
constrained geometrically, and has a higher average vi- 
brational frequency. 

Let us now come to E' activation. In the present pic- 
ture, E' is created transforming undistorted V Q into dis- 
torted Vq via two routes. The first proceeds on path A 



in Fig. with two successive one-photon ionizations of V Q 
into Vj 2 , followed by non-radiative decay into E'. This 
path is efficient since the +1 undistorted state is kept 
populated by sustained illumination (a much less effi- 
cient double-photon excitation of Vq into Vq 2 may also 
occur) . The excitation energies for path A are both near 
4-4.5 eV if the ionized electron is transfered to the Fermi 
level, i.e. to the Ep-pinning impurity; if it is promoted 
to the conduction band, the energies are instead about 
7 eV. Both processes are possible with X or 7 radiation, 
whose energy vastly exceeds that needed in the transi- 
tion. In UV irradiation, the center is often activated by 
pumping at 5 eV, and clearly only Fermi-level capture 
matches this figure. (There are, however, qualitative in- 
dications that the 5 eV excitation may activate E' via 
alternate routes involving other pre-existing defects.) 

The second excitation route, path B, involves an op- 
tical excitation of V Q into undistorted Vq , and a ther- 
mal excitation of the latter into the puckered state. The 
energy difference (0.3 eV) between the two Vq states 
implies that only a fraction of 1CP 5 of the vacancies 
gets promoted into the distorted state in equilibrium at 
room temperature, on sustained illumination. Therefore, 
though admissible, this path is preempted by path A. 
With the level ordering of Ref. ||] , also path B competes 
with path A. (Our ordering, however, matches better the 
thermal behavior, as discussed above.) 

Other calculated observables of our E' model are con- 
sistent with previous studies M-pLfUj. For instance, the 
optical absorption of the neutral undistorted state at 
6.9 eV correlates well with the 7.6 eV absorption band 
usually attributed to the neutral vacancy 0. For E' 
(metastable puckered Vq) we find an absorption at 4.7 
eV (defect-to-conduction promotion) followed by slow 
non-radiative decay back into E'. Since it is not followed 
by any emission or E' deactivation, this absorption cor- 
relates with the broad 5.8-eV band typical of E' ||il| , |i4| , 
which exhibits the same behavior in experiment. 

In summary, we conclusively put on firm ground the 
identification of the singly positive O vacancy in Si02 
with the E' center proving its thermodynamical stabil- 
ity via first principles calculations. Our picture provides 



naturally activation and deactivation mechanisms, and 
other optical signatures, in agreement with known exper- 
imental observations. In addition, our picture naturally 
explains the moderate p conditions produced by irradia- 
tion in Si02- 
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